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Zheng [Phys. Rev. E 61, 153 (2000)] claims that phase ordering dynamics in the microcanonical 
model displays unusual scaling laws. We show here, performing more careful numerical inves- 
tigations, that Zheng only observed transient dynamics mostly due to the corrections to scaling 
introduced by lattice effects, and that Ising-like (model A) phase ordering actually takes place at 
late times. Moreover, we argue that energy conservation manifests itself in different corrections to 
scaling. 



The problem of the dynamical foundations of statis- 
tical mechanics has received renewed attention recently, 
when, in the spirit of the famous work of Fermi, Pasta, 
and Ulam [1] , researchers directly studied the evolution of 
isolated, many-degrees-of-freedom Hamiltonian systems 
with the aim of relating their microscopic, deterministic, 
chaotic motion to their macroscopic statistical properties 
[2]. In this context, the two-dimensional lattice 0'* model 
is of special interest because it is known to exhibit, within 
the canonical ensemble, a second-order phase transition 
in the Ising universality class. Recent work [3,4] aimed, in 
particular, at studying the corresponding behavior in the 
isolated, microcanonical case, whose equations of motion 
can be written: 



where the sum is over the four nearest neighbors of site 
« on a square lattice. 

In the same spirit, Zheng [5] has considered the a 
priori simpler problem of the phase-ordering process 
which takes place when the microcanonical 0^ lattice 
is suddenly "quenched" below the critical point. Uni- 
versality of domain growth laws is nowadays a fairly 
well-established topic [6]. It is well documented, even 
for less traditional systems such as deterministic, possi- 
bly chaotic, spatially-extended dynamical systems [7-10]. 
For a scalar order parameter, two main universality 
classes can be distinguished depending on whether it is 
locally conserved or not. The Ising model is prototypical 
of the non-conserved case, and so should be the lattice (/)^ 
model, at least in the usual canonical ensemble point of 
view. However, and this was the interesting point raised 
by Zheng, the presence of the energy conservation in the 
microcanonical case might have an influence on the dy- 
namical scaling laws associated with phase ordering. In 
this sense, the question is whether the phase ordering of 
model C [11] is in the same universality class as model A. 

Using numerical simulations in the so-called "early- 
time" regime, Zheng confirmed that the usual dynam- 
ical scaling laws seem to hold, but with exponents at 
odds with those both the non-conserved order parame- 
ter (NCOP) and the conserved order parameter (COP) 



class [5]. In particular, he found z — 2.6(1) ( 1/z is the 
exponent governing the algebraic growth of L, the typi- 
cal size of domains), between its NCOP and COP values 
(^^NCOP = 2 and zcop = 3). Zheng "explains" this sur- 
prising result by the influence on scaling of "the fixed 
point corresponding to the minimum energy of random 
initial states" (present when considering domain growth) . 

Here we show that a more careful numerical investiga- 
tion actually leads to conclude that the microcanonical 
lattice (j)'^ model shows normal, NCOP phase-ordering. 
We argue that energy-conservation, but also lattice ef- 
fects, have sub-leading influences on scaling. We dis- 
cuss in particular the effect of the increase of the "bulk 
temperature" due to the progressive disappearance of in- 
terfaces during the growth transient. We attribute the 
erroneous results of Zheng to the danger of using "early- 
time" methods and naive logarithmic plots in problems 
with large microscopic times and/or corrections to scal- 
ing. 



I. NORMAL SCALING AT LATE TIMES 



Zheng conducted different numerical experiments 
which led to estimates of exponents z, X/z (where A is 
the so-called Fisher-Huse exponent), and 9 = (2 — X)/z. 
From all direct measurements of z he concluded that 
z = 2.6(1). Using this value, he found an estimate of A 
in agreement with the NCOP value (A = 1.22(5) whereas 
Ancop = 5/4)- Thus, the only strong departure from the 
NCOP values is for exponent z. Therefore, in the follow- 
ing, we focus on growth law of L in large systems at late 
times (as opposed to the early-time approach favored by 
Zheng) . 

In order to reach late times satisfactorily, we need a 
better control of the conservation of energy than with 
the simple second-order scheme used by Zheng. The fol- 
lowing results were obtained with a third-order bilateral 
symplectic algorithm [12] with a timestep dt = 0.025. 
The conservation of initial energy is better than 10^^ in 
relative value in all runs presented. To investigate phase- 
ordering, we use the same initial conditions as Zheng 
((/) = ±go where the sign is random and qo calculated 
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to yield the desired energy density e). We present results 
for two sets of parameter values, {rn^^g) = (6, 1.8) (set 
A, used by Zheng), and {m^.g) = (2,0.6) (set B, used in 
[3,4]). The initial energy density (eo = 27 for set A, and 
eo = 10.0001 for set B), was chosen very close to its mini- 
mum value allowed by the random-sign initial conditions 
(emin = 80/3 for set A, and emin = 10 for set B). This 
ensures that "thermal" fluctuations are minimized, since 
the energy density is then as far as possible from the crit- 
ical energy density Cc (cc ~ 35 for set A, and Cc « 21 for 
set B). 
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FIG. 1. Growth of typical scale L(t). All results presented 
here were obtained from single runs of lattices of 8192 x 8192 
sites with periodic boundary conditions. (a,b): L vs for 
parameter set A and B, the dotted lines arc the fits discussed 
in section II C; (c,d): same but in logarithmic scales (dashed 
lines for "expected" behavior z = 2). Insets: local slopes cal- 
culated by a running average of the derivated signal over a 
window of fixed size in the corresponding variable (i.e. \/t or 
logjQ t here). For (a,b) window size A(Vt) = 0.5; For (c,d) 
window size A(logio = 0.05. 



ing holds by observing the collapse of C(r, t) = f{r/L{t)) 
curves at different times after some transient (not shown) 
[19]. 

For both sets of parameters, the expected NCOP law 
(L ~ Vt, i.e. z = 2) is reached at late times, but a rather 
long transient is observed, especially for set A (Fig. lab). 
The same data plotted in logarithmic scales is thus mis- 
leading. If the data for set B reach the "normal" scaling 
(see Fig. Id and its inset) at late times, the correspond- 
ing plot for parameter set A (Fig. Ic) seems to indicate 
a value of 1 /z between i and ^ (typical of the value es- 
timated by Zheng) if one ignores the systematic trend 
upward of the local exponent (see inset of Fig. Ic) [20] . 
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FIG. 2. Decay of autocorrelation function A{t,to) from sin- 
gle runs of lattices of 4096 x 4096 as in Fig. 1 (to = 5). (a,b): 
Avst for parameter sets A and B (logarithmic scales, dashed 
lines: slope \/z = 5/8, insets: local slopes); (c,d): A vs L 
(logarithmic scales) for parameter sets A and B (logarithmic 
scales, dashed lines: slope A = 5/4, the quality of the data is 
not high enough for local slopes). 



A. Growth of typical domain size 

The typical domain size L was determined by the mid- 
height value of C{r,t), the normalized two-point corre- 
lation function calculated for simplicity along the prin- 
cipal axes of the lattice using the continuous field (p 
or the reduced "spin" variable a = sign((/)): C{k,t) = 



2N Si 



+ k ] 



(7, 



i j + k) 



(No significant difference 



was found between the two cases, and only results using a 
are shown below.) We first checked that dynamical scal- 



B. Decay of autocorrelation function 

Following Zheng, wc also measured the decay of the au- 
tocorrelation function A{t,to) = {a{t)a{to)), but, again, 
for rather large systems and late times, choosing, in par- 
ticular, initial reference times to larger than the "micro- 
scopic" transient time. In the usual NCOP phase order- 
ing framework, we expect 



A{t,to) 



L{to) 
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where the Fisher-Husc exponent A = 5/4. Again, plot- 
ting log A vs log t for parameter set A may yields an "ef- 
fective" exponent X/z smaller than its NCOP value (and 
close to the value given by Zheng), but a closer look re- 
veals a systematic increase of the instantaneous exponent 
(Fig. 2a and inset). However, plotting A vs L, the ex- 
pected scaling is observed (Fig. 2c). Considering now 
parameter set B, the expected scaling laws are observed 
rather easily (Fig. 2b, d). This confirms further that for 
the parameter values chosen by Zheng, the onset of the 
asymptotic scaling regime is delayed. 



(a) (b) 




FIG. 3. Early-time growth of the squared magnetization 
for an ensemble of 500 symmetric {M = 0) samples in a lat- 
tice of 512 X 512 sites. (a,b): (M^) vs t (logarithmic scales) for 
parameter sets A and B; dashed lines: slope 1. (c,d): (M^) 
vs t (linear scales) for parameter sets A and B. 



C. Early-time growth of squared magnetization 

For the sake of completeness, and in order to probe 
the validity of the early-time scaling approach taken 
by Zheng, we also performed numerical simulations to 
measure the short-time growth of the squared magneti- 
zation (M being defined as the spatial average of 
either or ct, = sign((^j)) for samples with zero ini- 
tial magnetization. In this case, for NCOP scaling, we 
expect (M^) ^ f^^^ = t. Our data on logarithmic 
scales barely reaches the expected behavior at f ~ 10'^ 
(Fig. 3a,b). Note that the corrections have different sign 
for the two parameter sets. On linear scales, however, our 
data reveals the expected proportionality of (M^) and t 
(Fig. 3c,d). 



II. CORRECTIONS TO SCALING: ENERGY 
CONSERVATION AND LATTICE EFFECTS 



The above results show that domain growth in the mi- 
crocanonical model eventually falls into the NCOP 
universality class {z = 2, A = 5/4) after some possibly 
long transient behavior. In this section, we suggest that 
there are two main factors at the origin of these correc- 
tions to scaling: space discretization and energy conser- 
vation. 



A. Lattice effects 



For the sets of parameters studied by Zheng (notably 

parameter set A), domain growth initially appears to be 
slower than the expected NCOP law (the effective value 
of z measured at short times is larger than two). 

This is similar to earlier observations on coupled map 
lattices, both for the NCOP and COP cases. First mea- 
surements of domain growth seemed to indicate slower 
growth in those discrete-space, discrete-time chaotic 
models [7,9], but it was shown later that this non-trivial 
scaling is only transient and that, for late-enough times, 
normal scaling is recovered [8,10]. It was also argued 
that these long transients gradually disappear in the 
continuous-space limit which is well-defined in these sys- 
tems [13]. 

Here, it is easy to observe that, increasing g (and m?, 
since there is only one free parameter) , the slowing-down 
of domain growth due to lattice effects diminishes, as sug- 
gested by the difference between parameter sets A and 
B. For parameter set B, the transient actually presents 
faster growth than asymptotically: corrections to scaling 
are dominated by another phenomenon, stronger than 
lattice effects. 
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FIG. 4. Corrections to scaling due to energy conservation. 

Results for parameter set B, for which these corrections are 
dominant (see text), (a): the three components of the total 
energy during phase ordering, (b): same but from ordered 
initial conditions, (c): decay on the number of interfacial 
sites (defined as sites with a least one nearest-neighbor of op- 
posite sign); the number of interfacial sites measured during 
the evolution from an initially-homogeneous configuration has 
been substracted. (d): Convergence of the total kinetic en- 
ergy (temperature) towards its asymptotic value: 1/ Aitkin vs 
y/t during phase ordering. 



B. Energy fluxes 

The (conserved) total energy in the system can be par- 
titioned into three components, kinetic, potential, and 
interaction energy: 

El • 2 
2 

i 

i 

1 

i j = l 

As domain growth proceeds following random initial con- 
ditions, the interaction energy decreases with time as a 
large part of it is contained in the interfaces separating 
domains, the density of which decreases like 1/Vt. Po- 
tential decreases and kinetic energy increases (Fig. 4). 

Starting from "ordered" initial conditions (all sites in 
the same phase, e.g. Vi (/'i(O) > 0), the three components 



of the total energy are almost constant, except for a slight 
(quasi logarithmic in time) increase of interaction energy 
and decrease of potential energy due to the nucleation of 
"thermal" droplets and to relaxation further into poten- 
tial wells. This effect is indeed larger for parameter set 
A for which the minimal possible energy density emin is 
further away from the critical energy density, leading to 
stronger thermal fluctuations than for parameter set B. 



C. Corrections to effective temperature 

Kinetic energy can be interpreted as a temperature in 
the microcanonical context [3]. We can thus see the ob- 
served increase of kinetic energy as an increase of the tem- 
perature of the system. For the (canonical) Ising model, 
it is well known [14] that the prefactor of the growth law 
of L(t) decreases to zero as the temperature approaches 
its critical value. We suggest to see the growth law ob- 
served here for the ip'^ model as including a "temperature- 
dependent" prefactor: 



L{t) ~ K{E^i,,) Vi 



(1) 



Quantitatively (Fig. 4), the kinetic energy seems to reach 
its asymptotic value like: 

AEkin = Skin(i) - i^ki„(oo) ~ (2) 
Assuming its analyticity, we can write the prefactor K: 

K{EMn) ^Koo+K[ ASkin + (A^kin)' + . . . (3) 

where Koo is the asymptotic {t — > oo) prefactor of the 
domain growth law. Injecting (2) and (3) into (1), we fi- 
nally expect the following Ansatz to hold for the domain 
growth law: 



El 



(4) 



Equation (4) provides excellent fits to our data for the 
growth of L{t). We find K^o ^ 0.2, Ki ~ 2.0, and 
K2 — —3.5 for parameter set A (with ~ 7.8 and cor- 
relation coefficient 0.99996) , and Koc — 2.1, Ki ~ —5.6, 
and K2 — 5.7 (with = 10^ and correlation coefficient 
0.9997) for parameter set B. 

The corrective terms have opposite signs in both cases: 
this indicates that Eq. (4) is only relevant for parameter 
set B, because, in analogy with the Ising model, only 
negative values of Ki are allowed (the prefactor of the 
growth law decreases with increasing temperature). A 
first conclusion is thus that the main corrections to scal- 
ing for parameter set A (typical of those used by Zheng) 
are due to lattice effects and are not consequences of the 
conservation of global energy. We note that, similarly, a 
fit of domain growth in coupled map lattices also yields 
positive values of Ki , indicating that this sign is a signa- 
ture of lattice effects [8] . 
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On the other hand, the above framework does provide 
the relevant explanation for the corrections to scaling ob- 
served for parameter set B, which can thus be traced 
back to the fluxes between the various components of 
the energy induced by the decrease of interfaces between 
domains as phase-ordering proceeds. 

III. CONCLUSION 

In this Comment, we have shown that Zheng reached 

erroneous conclusions when studying phase ordering in 
the microcanonical lattice 0^ model. This system, like 
other chaotic, deterministic, dynamical systems present- 
ing phase ordering, does display the expected domain 
growth scaling laws, i.e. those of the non-conserved or- 
der parameter case (model A, z = 2, A = 5/4). We have 
shown further that the main influence of the conservation 
of energy is to introduce corrections to scaling, but that 
the long transients which plagued Zheng's approach are 
due to lattice effects. 

Zheng offered, as an explanation for his non- 
conventional results, that the system considered falls into 
the class of model C, where a conserved density is coupled 
to the order parameter [11]. Studies of phase-ordering 
in model C [16,17,15] show that model B-like, but also 
model A-like behavior can be observed. In this con- 
text, the microcanonical lattice can be considered as 
a model C system quenched into its "bistable region" 
(where 2; = 2 is observed [18,15]). 

We thank A. Lemaitre and A. S. Somoza for interesting 
discussions and L. Casetti for providing a sample version 
of his program. 
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